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Abstract 



Recent concept of integrable nonholonomic deformation found for the KdV equation is 
, extended to the mKdV equation and generahzed to the AKNS system. For the deformed 

I mKdV equation we find a matrix Lax pair, a novel two-fold integrable hierarchy and 

■ exact N-soliton solutions exhibiting unusual accelerating motion. We show that both 

. the deformed KdV and mKdV systems possess infinitely many generalized symmetries, 

O I conserved quantities and a recursion operator. 
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1 Introduction 



Recently discovered integrable 6-th order Korteweg de Vries (KdV) equation [T] was shown 
to represent a nonholonomic deformation (NHD) of the well known KdV equation preserving 
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its integrability and exhibiting an integrable hierarchy In a subsequent development a 
matrix Lax pair, the N-sohton solution through inverse scattering transform (1ST) method 
and an intriguing two-fold integrable hierarchy were found for this particular system by one 
of the authors |3] . The NHD of the KdV (dKdV) was shown recently to be a certain form of 
self-consistent source equation (SCSE) allowing particular exact solutions |4]. 

The NHD for such field theoretical models is given by a constraint in the form of nonlinear 
differential equation involving only x -derivatives on a single perturbing function, which is 
deforming the original integrable equation. This type of integrable deformation is a relatively 
recent discovery, which allows also an integrable hierarchy of higher order deformations. On the 
other hand the construction of SCSE is a well known concept [8j, which is represented usually 
by a coupled system consisting of the original integrable equation with an additional term 
made from the eigenfunctions and 2N eigenvalue equations of the Lax operator with explicit 
dependence on N discrete eigenvalues. In general the eigenfunctions and eigenvalues are complex 
and one can also add another 2N number of evolution equations for the eigenfunctions. In a 
recent development of such SCSE its Rosochatius deformation is proposed [3]. At certain limit 
or at particular cases the contact between the lowest order NHD and the SCSE can be found, 
though the integrable hierarchy of higher order deformations seems to be possible only for the 
NHD. 

In a slightly generalized form this deformed KdV equation can be given as 

gxxx + 4:ugx + 2ux{g + c{t)) = 0, (2) 

where subscript denote partial derivatives. Though the time-dependent arbitrary function c{t) 
can be absorbed in equation ([2]) by redefining perturbing function g{x,t) , we keep it in the 
explicit form for our convenience. Recall that, the KdV and the modified KdV (mKdV) are 
intimately related partner systems [5|. Therefore it would be natural to expect that the concept 
of NHD found for the KdV should also be extendible to the mKdV equation. We show here that 
this expectation is indeed true by constructing explicitly a novel integrable NHD of the mKdV 
(dmKdV) equation, which can yield a new integrable 4th order potential mKdV equation. 
We discover for this integrable deformed mKdV equation with a nonholonomic constraint a 
matrix Lax pair, exact N-soliton solutions and a novel two-fold integrable hierarchy, similar 
to the result of the dKdV [3]. The solitons, found for both the basic field and the perturbing 
function of the dmKdV, show unusual accelerating (or decelerating) motion. We show that 
both the deformed KdV and mKdV systems possess infinitely many generalized symmetries 
and conserved quantities, and a recursion operator which have been studied until now only for 
their undeformed counterparts. The Lie symmetry analysis we perform here for the deformed 
equations leads also to explicit construction for the hierarchies of the generalized symmetries 
for both the basic field and the perturbing functions. Our analysis of the continuity equation 
for the dKdV and dmKdV reveals that the conserved densities for these deformed systems 
remain the same as in their original undeformed cases, whereas the current densities (fluxes) 
explicitly contain the deforming functions. This shows the intriguing fact that the nonholonomic 
deformations can appear only at the equation level, while the conserved integrals of motion 
remain the same under NHD. 
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Finally we unify dKdV and dmKdV to discover an integrable nonholonomic deformation 
for the more general AKNS [5j system. 

The plan of the paper is as follows. In section 2 the new NHD of the mKdV equation 
and the related integrability structures such as the matrix Lax pair and the exact A^-soliton 
solution are presented. Tt is also shown that the dmKdV admits infinitely many higher order or 
generalized symmetries, conserved quantities and a recursion operator. In section 3 a similar 
analysis has been carried out for the dKdV equation. Section 4 generalizes the nonholonomic 
deformation to the AKNS system, revealing a novel two-fold integrable hierarchy for all its 
members and consequently for the dmKdV and the dKdV systems. In section 5 we give a brief 
summary of our results and the concluding remarks. 



2 Nonholonomic deformation of the mKdV equation 

Integrable equations with NHD should be driven by an additional perturbative or deforming 
function, which in turn would be subjected to a differential constraint of nonholonomic nature. 
Therefore, analogous to the deformed KdV equation ([T]l2l) we propose a deformed modified 
KdV equation with nonholonomic constraint as 

vt - fxxx - Gv^Vx = w{t, x), (3) 
w^-2v{c'^{t) -w^)^ = 0. (4) 

Note that the arbitrary function c{t) can be removed again from the equation (jlj) by rescaling 
^^^^^ w{x, t) . We however keep c{t) in the explicit form for the later convenience. At 
c(t) = we should have the deforming function w{x,t) = 0, when the dmKdV (l3}|4]) would 
reduce to the standard mKdV equation. For establishing the integrability of (JSHll) we derive 
the associated pair of matrix Lax operators in the form 

[/(A) = U"''"^''iX), y(A) = y'"'='^"(A) + V^'^iX), (5) 
where U"^^'^'"{X), V"^'^'^^{X) are the well known Lax pair for the standard mKdV system [5]: 
Umkdv^^-^ = iXa^ + iva\ 

ymkdv^^-^ = 2iX{v^ - 2X)a^ - 2iXv^a^ + i(v^^ - AX'^v + 2v^)a\ (6) 
while the additional term 

V^'^iX) = -\{ba^ - wa^), where b = {c\t) - (7) 
2X 

is responsible for the deformation of the mKdV equation. Here a"", a = 1,2,3 are standard 
Pauli matrices 

^ = ^ 1 J ' ^ = ^ z J ' ^ = U -1 

and = |((T^ ±ia^). For confirming the Lax integrability of the NHD of the mKdV equation 
we show that it is derivable from the flatness condition Ut — + [U,V] = of the 
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Lax pair In this process one finds that, the expressions in all positive powers of A",n = 
1,2,3 vanishes trivially, while that for n = yields the deformed niKdV equation . The 
coefficients with the power n = —1 gives in turn the constraint equations 

bx = —2vw, Wx = 2vb (9) 

from which eliminating function b by using b = (c^(t) — w"^)^ we arrive at the constraint (jl]) 
as required. 

Note that one can also find a single higher order nonlinear equation by eliminating further 
the deforming function w{x, t) from the set of equations ([3lll]). By introducing a potential field 
9x = 2v we can rewrite Eq. (jlj) in the form 6^ = — — — r, giving an easy solution w = c{t) sin 6 , 

(c2--u;2)3 

which can be checked by direct substitution, since Wx = c6x cos 9 while (c^ — w'^)2 = ccos9 . 
Inserting these expressions of w , v and hence those of Vt, Vxxx, v'^Vx through 9 in Eq. one 
can derive a new integrable 4-th order potential mKdV equation in 9 given by 

{9t~9xxx-\9l)x = 2cit)sm9. (10) 

Interestingly, keeping only the term 9xt in the LHS, while grouping the other terms as 

f{9) = {9xxx + lel)x, (11) 

one can rewrite Eq. ffTOj) in the form of a perturbed sine-Gordon (SG) equation in the light-cone 
coordinates as 

^,, = 2c(t)sin^ + /(^), (12) 

with a time-dependent mass parameter = 2c{t) , perturbed by a function fllip . Significantly 
unlike most of the perturbed SG equation fll2tilip . derived here, is integrable allowing exact 
N-soliton solution. Thus the integrable higher order dmKdV equation ffTOl) is equivalent to the 
NHD of the dmKdV equation ([3lll]), which in turn is equivalent to the set of equations ([3]) with 
©• 

2.1 Exact soliton solutions 

As is well known any perturbation usually spoils the integrability of a nonlinear system and 
hence forbids general analytic solutions. However the dmKdV we have constructed retains its 
complete integrability in spite of the perturbation with nonholonomic constraint. Therefore ex- 
ploiting this integrability property we intend to derive exact N-soliton solutions for the dmKdV 
equation or equivalently for the novel 4-th order potential mKdV equation flTUl) . through 
the 1ST method. The procedure follows that for the standard mKdV equation [6] in its initial 
steps, while in the final step the effect of deformation should be incorporated. It is intriguing 
to note from ([5]), that for the dmKdV equation, only the time Lax operator V is deformed, 
while the space part U is kept unchanged. This reveals an important connection between the 
deformed time evolution of the Jost solution and the NHD of the nonlinear equation. 
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Recall that for extracting the exact soliton solution in the 1ST the space-Lax operator U (A) 
describing the scattering of the Jost functions, plays the key role. Only at the final stage we 
need to fix the time evolution of the solitons through the time- dependence of the spectral data, 
which in turn is determined from the asymptotic value of the time-Lax operator l^(A)||a;|^oo • 
Therefore following [6] we can derive the N-soliton solution for our deformed mKdV equation 
as 

cf 1 
vix) = [—^iTidetAix)]^ , (13) 

where the matrix function A{x) is expressed through its elements as 
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Anm = J2 f-lflm, whcrC fnm = Kn. + ^^e-(''"+'^™)^ (14) 

Here parameters K„,n = 1,2,. ..,A^, denote the time- independent zeros of the scattering 
matrix element: a(A = A„) = 0, along the imaginary axis: A„ = and /3n(^) = f'(A = A„) 
are the time-dependent spectral data to be determined from K(A) = V'^^'^'"[X) -|- V^^^ [X) , at 
X — > ±00 . 

Note that due to the boundary condition (BC) ti— s>0,iy— >0 at x— ±00 , the asymptotic 
value of ([6]): V'^^'^'"[X) — >■ —4iA^ cr^ corresponds to the undeformed part, while the BC h — >■ c(t) 
at X ^ ±cxD , determines the crucial effect of deformation y^'^'^(A) |A^^c(t)a^ . As a result 
we obtain 

Ut) = /?„(0)e(^''"*+^), ctit) = c(t), (15) 
yielding finally the {x,t) dependent exact soliton solution from (IHj). 

To see the effect of deformation on the dynamics of the solitons more closely we construct 
1-soliton solution for the dmKDV equation as reduced from (I13H15I) at = 1 : 

v{x, t) = -ysech^, ^ = k{x — vt) + (f), v = vq + Vd (16) 

where the phase (p is an arbitrary constant, Vq = is the usual constant velocity of the 
mKdV soliton, while Vd = is the unusual time-dependent part of the velocity, induced by 
the deformation. Note that the time- dependent asymptotic value of the deformation acts here 
like a forcing term sitting at the space boundaries, which for c{t) = c^t with cq > forces the 
soliton to accelerate, while with cq < makes it decelerate (see Fig 1). 

It is important to note, that the evolution of the basic field v{x,t) comes here from two 
sources: dtV = dt^v + dt^v , induced by two different 'times' to and td- dt^v = dtv\c{t)=o is 
the evolution due to the unperturbed time to , caused by the standard dispersive and nonlinear 
terms in ([3]), while dt^v = c{t)dc{t)V is the evolution due to the deformed time , caused by 
the perturbation linked as 

w = dt^v, b = c{t) + / vdt^vdx. (17) 



Therefore using (ITTI) and dc{t)f^nif) = ^n^n{t) we can find from fl 13111 51) the N-soliton solution 
corresponding to w,b, which for A = 1 derived from f|T6|) gives also localized accelerating 
solution 

w{x,t) = c(t)/«sech^tanh(^ and b{x,t) = c(t)(l — «;^sech^^). (18) 
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This shows that the perturbing function, taking itself the sohtonic form drives the field soliton 
to have an accelerated motion, while in turn the basic field solution self-consistently determines 
the solitonic form of the perturbing function, sustaining thus the integrability of the system. 



2.2 Generalized symmetries 

It is well known that the Lie symmetry analysis plays an effective role to study the integra- 
bility properties of nonlinear evolution equations in (1+1) dimensions such as the existence of 
infinitely many generalised symmetries, conserved quantities and a recursion operator ([9], [TO] . 
|llj . [T2] . |13j). We show here that similar analysis can be performed with equal success for 
our deformed mKdV equation 



Note that the dmKdV equation is invariant under the scaling or dilatation symmetry 

{t, X, V, w, c(t)) {s~'^t, s~^x, s'^v, s'^w, s^c{t)), 

where s is an arbitrary parameter which suggests that v corresponds to one derivative with re- 
spect to scaling x , w corresponds to three derivatives with respect to x and c{t) corresponds 
to four derivatives with respect to x . We would like to mention that Hereman and his collabo- 
rators have developed an algorithm to derive generalised symmetries, conserved quantities and 
recursion operators for nonlinear partial differential and differential-difference equations [H]. 
Hereman's algorithm is based basically on the concept of weights and ranks. The weight W of 
a variable is defined as the exponent in the scaling parameter s which multiplies the variable. 
Weights of the dependent variables are non-negative and rational. An expression is said to be 
uniform in rank if all its terms have the same rank. Setting W{D) = W{£) = 1 , we see that 
W{v) = 1, W{w) = 4, W{b) = 4, W{c{t)) = 4 and W{-§-^) = 3 and hence eqn. ([3D is of rank 4 
and (jl]) is of rank 5 . This property is called the uniformity in rank. The rank of a monomial 
is defined as the total weight of the monomial, again in terms of derivatives with respect to x . 

Now, assume that the deformed mKdV equation ([3}|1]) is invariant under one parameter 
nonpoint continuous transformations 

t* = t, x* = x, V* = v + eGi + 0{e^), w* = w + eHi + 0{t^), (19) 

where 

Gi Gi(^Vj W, Vxi Vxxi '^xxxi '^4x5 Wxxj ' ' ') 

Hi Hi (f, If , Vx , Vxx ) "^xxx 1 V^x ; • • ; '^x i '^xx ; ' ' ') 

provided v{x,t) and w{x,t) satisfy equation ([3]- H]). Consequently we obtain the following 
invariant equations 

DGi D'^Gi ^DGi ^ DHi , , 

Dt Dx^ Dx Dx ' ^ ^ 

—i = 2G.Vic{ty - w^) - 2vw^ , (21) 
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Gi = Vx, Hi = w, (22) 

is a trivial generahzed symmetry with rank (2,4). This suggests that the next generahzed 
symmetry G2 and H2 of dmKdV must have rank (4,6). With this in mind we first form 
monomials in v and w of rank (4,6), Thus the most general form of G2 and H2 will be 

G2 = aiv'^Vx + a2V3x + asw, H2 = biv'^w + b2W2x, (23) 

where ai,bi,i = 1,2,3 are arbitrary constants to be determined. We now substitute G2 and 
H2 in the invariant equation ( !20ll2Tl) .with i = 2 and using (ISM]) we find that the consistency 
condition holds only for the following parametric restrictions: 

fli = 6, 02 = 1, ^3 = 0, 61 = 6, ^2 = 1 
and so the generalized symmetry with rank (4, 6) becomes 

G2 = Qv'^Vx + Vsx, H2 = Qv'^W + W2x- (24) 

Proceeding as above, for z = 3 , we find that the invariant equations (12011211) satisfy only if 

G3 = V^x + + 40vVxV2x + I0v\3x + SOv\x, (25) 
H3 = Wix + Wv^W + 20vV2xW + 20vVxWx + l0v'^W2x + 30f ^1(7, (26) 

which is a next nontrivial generalized symmetries with rank( 6,8). In a similar manner we 
can derive infinitely many generalized symmetries {(G4, H4), (G5, H^), ■ ■ •} for ([3111]) with rank 
{(8, 10), (10, 12), ■ ■ ■} . We have also checked that the commutator 

[G„Gi+,] = 0,[H„H,+,]=0, Vz, 

indicating that dmKdV eqn. ([3111]) admits an infinitely many generalized symmetries which 
commute. 



2.3 Recursion operator 

In this section we derive a Recursion operator, an important property of the integrable systems, 
for our deformed mKdV eqn. ([Hill]), which is usually possible to obtain only for the unperturbed 
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integrable systems. An operator valued function TZ is said to be a recursion operator of a scalar 
nonlinear partial differential equation with two independent variables if it satisfies 

G = nG, 

where G and G are successive generalized symmetries. For the above equation can be 
written as 

Gk+i 
Hk+i 

where {Gk^H^) and {Gk+i, Hk+i) are successive generalized symmetries and Rij,i,j = 1,2 are 
functions of dependent variable and their differential and integral operators. The construction 
of the recursion operator TZ for the dmKdV equation is as follows: For k = 2 equation (1271) 
becomes 

_ Hs 

where (G2,-f^2) and (Gsji^s) are the generalized symmetries of ranks (4, 6) and (6,8) respec- 
tively. The ranks of -R12, i?2i and R22 can be determined from the following relations 

rankG^ = rankRu + rankG2 = rankRi2 + rankH2, (29) 

rankH^ = rankR2i + rankG2 = rankR22 + rankH2, (30) 

Equations fl29tl30p show that the ranks of Ru, R12, R21, R22 respectively are 2, 0, 4, 2 and so we 
consider the entries of TZ written in terms of differential and integral operators of the dependent 
variables having the form 

-Rll = Co(9^ + civ'^ + C2V^d'^v, R12 = 0, R22 = god'^ + 9iv'^ + 92vd~^v^, 

R21 = fowd-'v + f,vd-'w + f2v' + fsw + Uv,,d-' + f,v^d\ (31) 

where Ci, gi,i = 0,1,2 and fj,j = 0,1,2,3,4,5 are constants to be determined. Substituting 
the above in fl28p we find that it is satisfied identically only if 

Co = 1, ci = 4, C2 = 4:Jo = 4, /i = 4, /2 = /a = /4 = /5 = 0, ^0 = 1, ^1 = 4, 92 = -4 
and therefore the recursion operator for the dmKdV equation becomes 

^_ f d'^ + 4v^ + Av^d-^v 



TZ 



Gk 
Hk 



RU R12 
R21 R22 



Gk 
Hk 



(27) 



Rn R12 
R21 R22 



G2 
H2 



(28) 




2.4 Conserved quantities 



A local conservation law of a partial differential equation with two independent variables (x, t) 
is defined by 
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which is satisfied on all solutions. The function p{x, t) is usually called local conserved density 
and J(x, t) is the associated fiux also known as current density. We show that the dmKdV 
([SlH]) or eqn.(3) and (8) admit infinitely many polynomial conserved quantities. From ([3]), we 
find directly that 



p[l] = V, J[l] = -{v,, + 2v'+ / dxw) (34) 



is a trivial conserved quantity with rank (1,3). Recall that ([3]) and (8) is invariant under 
scaling symmetry {t,x,v,w,b,c) {s~^t, s~^x, s^v, s'^w, 3%, s'^c), where s is an arbitrary 
parameter. To derive a conserved quantity with rank (2,4) , as before, we form monomials of 
v{x,t) and w{x,t) which gives the list £i = {v'^} . Thus the conserved density of rank 2 will 
be p{x,t) = f^. As a result we obtain 

p[2] = ^;^ J[2] = -2vv^^ - 3v^ + vl + b. (35) 

Proceeding as above we find the next two conserved quantities as 

p[3] =v^- vl, J[3] = -Av'' - Avh,, + 2v,v^^^ - vl, + Uv^l + 2v\ (36) 

p[4] = v^-5v'vl + ^vl, (37) 
—9 1 

^[4] = —V^ - Qv^V^cx + IQv^VxVxxx - VxxVix + ^^Lx " S^^^L 

+ 4:5v'^vl-10vvlv:,^-^v^ + 3v%-vlb-2vv,b^. (38) 

with ranks (4,6) and (6,8) respectively. In a similar manner we can derive an infinitely many 
conserved quantities {p{x,t), J{x,t)) for dmKdV with ranks {(8, 10), (10, 12) . . .} which in- 
volve lengthy expressions and so the details are omitted here. 

It is important to notice that the integrals of motion describing infinite number of in- 
tegrated conserved quantities Cn, n = 1,2, . . ., which should be commutative as a necessary 
criterion of Liouville integrability can be given by c„ = J p[n]dx . It is clearly seen from the con- 
tinuity equation (1531) that due to vanishing of the fields along with their derivatives at the space- 
infinities we naturally obtain dfCn = . Therefore from the expressions of p[n], n = 1, 2, . . . , 
derived above we get the conserved quantities for the dmKdV equation as 

(39) 

with C3 = Hmkdv being the Hamiltonian of the mKdV equation. We see therefore that the 
conserved quantities including the Hamiltonian remains the same for both the deformed and 
undeformed mKdV systems, though the corresponding equations for the dmKdV with addi- 
tional perturbing function and nonholonomic constraint on it are surely different. Note also that 
this effect of deformation changes the structure of the local current densities J[n], n = 1, 2, . . . , 
which contain the deforming functions but not the densities p[rt\ , which generate the 

conserved quantities. 
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3 Nonholonomic deformation of the KdV equation and 
its integrability aspects 



Nonholonomic constraints on field models have received increasing attention over recent years 
, while a significant breakthrough is made recently by discovering an integrable nonholonomic 
deformation for the KdV equation [H [2|, [3l H] . It has been established that this deformed KdV 
equation like the undeformed standard KdV admits Lax pair, N-soliton solutions and a two-fold 
integrable hierarchy [3]. The deformed KdV can be considered also as a source equation [1], 
which is however different from and simpler than the well known source KdV equation [8] . The 
novelty of this source is that it can be deformed recursively by going to the next order in its 
integrable hierarchy with higher order deformations [31 [7]. 

However many other fundamental and important properties of the integrable systems, 
such as the existence of an infinite number of generalized symmetries, conserved quantities and 
a recursion operator, which have been well established for the standard KdV equation have 
not yet been studied for their nonholonomically deformed extension, investigation of which is 
therefore our aim here. 



3.1 Generalized symmetries 

The deformed KdV equation ([Tll2]) is obviously invariant under the dilatation symmetry 

{t, X, u, g, c(t)) {s~^t, s~^x, s'^u, s^g, s'^c{t)), (40) 

where s is an arbitrary parameter which suggests that u corresponds to two derivatives, 
while both g and c{t) correspond to four derivatives, with respect to x . Setting W{D) = 
W{£) = 1 , therefore one gets W{u) = 2, W{g) = 4, W{c{t)) = 4 and W{§-J = 3 and hence 
equations is of rank (5, 7) . This property is called the uniformity in rank. Now, assume 
that the deformed KdV equation is invariant under one parameter continuous nonpoint 
transformations 

u* = u + eG., + Oie'), (41) 
g* = g + eH, + 0{e^), (42) 

where 

C^i Gi{u^ g, U'xxj ^xxxy ^4x) ■•) 9x1 9xxi ' ' ') 
Hi Hi{u^ g, Ux, Uxxi Uxxxy U4.X1 ■•) gxi gxxi ' ' ■) 

provided u{x,t) and g{x,t) satisfy equations ([T]l2]). Consequently we obtain the following 
invariant equations 

DG, D'G, DG. DH, 

= T- + 6tt— h GuxGi + — - — , 43 

Dt Dx^ Dx Dx ^ ' 

^'^^ + + AG^gx + 2^{g + c(t)) + 2uxH, = 0. (44) 



Dx^ Dx Dx 
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From 



we derive the following generalized symmetries 

Gi = Ma;, Hi = G2 = Quux + u-ix, H2 = Qug^ + g^^, 



(45) 



which are trivial ones with ranks (3, 5) and (5, 7) , respectively. Proceeding as before we find 
that (P- [2]) admits an infinitely many generalized symmetries, where the first two nontrivial 
generalised symmetries are 



= 30u^g^ + IQu^gxx + Wu2xgx + ^OuQsx + 95x, 



and 



Hi = UOu^g^ + l^ulg^ + lAQuu.^g2x + l4Quu2xgx + 70u^g3x 
+4:2u2xg3x + 28u3xg2x + 2Su^g4x + 14:U4xgx + I4:ug5x + grx 

with ranks (7, 9) and (9, 11) , respectively. Note that the remaining higher order generalised 
symmetries involve lengthy expressions and hence not presented here. We find also that the 
commutator relations 

[G^G.+i] =0, [if„/f,+i] = 0,Vz 
hold. It is straightforward again to check that the derived generalised symmetries satisfy 



Hi+i 



n 



G, 
H 



Vz, 



with the recursion operator 

7^ = 



92 + 4u + 2uxd-^ 
2g^d-^ + 2d-^g^ d'^ + Au - 2d~^u,^ 



(46) 



3.2 Conserved quantities 

Recall that the dKdV equation ([T]l2]) is invariant under the dilatation symmetry [t, x, u, g) 
{s~^t, s~^x, s'^u, s^g). We can show again that the dKdV equation admits an infinitely many 
conserved quantities. From ([I]l2]), we find that 

p[l] = u, J[l] = -{u,, + 3u^ + g) (47) 

correspond to a trivial conserved quantity with rank (2,4). To derive the next conserved 
quantity with rank (4, 6) we form monomials of u and g as before, which give the list £1 = 
{u^, Uxx} ■ Thus the most general form of the conserved density of rank 4 would be p{x, t) = 
CiV? + C2Uxx, with ci , C2 are constants. It is straightforward to check that the next order 
conserved density and associated flux 

p[2] = u\ J[2] = (-4m^ - 2uu^x + ul + 2ug + gxx + 2uc{t)) (48) 
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satisfy fl33p . In a similar manner we obtain in the next higher order 

p[3] = (4m^ - 2ul), J[3] = {-18u^ + 24uul - 12u\^^ - 2ul^ + Au^u^^ 

- g4x + 4:U^g + 4:U^c(t) - 2uxxc{t) - 2u^^g - 6u^g^). (49) 

We wish to add that the next higher order conserved quantities involve lengthy expressions and 
so we refrain from presenting them. 

We can obviously construct the commuting family of integrated conserved quantities as 
above: c„ = / p[n]dx with dtCn = . Therefore from the corresponding p[n\,n = 1, 2, . . . , we 
get the conserved quantities for the dKdV equation as 

Ci = J u dx, C2 = j dx, ''3 ~ y ^^''-'^ ~ '^'^l)dx, ... (50) 

with C3 = Hkdv being the Hamiltonian of the KdV equation. Therefore we see again that for 
the deformed KdV system, the conserved quantities remain the same as in the undeformed 
case, though the corresponding equation gets deformed due to the nonholonomic constraint. 



4 Generalization to the nonholonomic deformation of 
the AKNS system 

In sect. 2 we have discovered a new integrable nonholonomic deformation of the mKdV system, 
extending the dKdV. Here we intend to show that these deformed systems can be unified and 
generalized further to an integrable nonholonomic deformation of the AKNS [5] system: 

qt - Qxxx - G{qr)qx = 9x, (51) 

n - Txxx - 6(gr)r^. = (52) 

For finding the integrable nonholonomic constraints on the deforming functions g, f we in- 
troduce the deforming matrix C^^^ = ga^ + fa~ + ha^ and similarly for G'-^-* , and de- 
note the AKNS fields in the matrix form ^Z'-*'-* = qa^ + ra^ . Integrability condition, i.e. 
the flatness condition of the associated deformed Lax pair (with the additional deformation 
V2'^^^\X) = ^{X^^G^^^ + A^^G*^^-*-* ), therefore leads to the nonholonomic constraint: 

Gi^) = z[[/(°),G«] +z[a3,G(2)], Gf = i[U^'\G^% (53) 

Higher order integrable deformations can be generated recursively by adding to V2 ^'^'^\X) more 
terms like ^\~^G^^\ j = 3, . . .n with arbitrary n . It would result to a new integrable hierarchy 
of nonholonomic deformations of the AKNS (dAKNS) system, given explicitly through the 
higher order constraints 

G« = z[^(°),G«]+z[a3,G(^)], ... , 
G("-i) = z[?7(°\G("-^)] + 2[a3,G(")], 
Qin) ^ z[[/(o),g(")]. (54) 
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For n = 1 this hierarchy clearly reduces to 

Gii) = 2[[/(°),G«], (55) 

while for n = 2 gives (!53|) . 

The integrability of the dAKNS system fl5m53p is guaranteed from its associated matrix 
Lax pair 

U{X) = iXa^ + iU'^^\ U^^'^ = qa+ + ra- , 

V{X) = y^^^^(A) + y^^^(A), (56) 

where 

V^™(A) = zUi^J - + 2a3A(-[/f + ^([/(°))2) - 4^[/Wa2 

+ 22(f/W)3-[t/W,[/f ], 

V^'^iX) = i(A-iG'W + A-2G(2)). (57) 

We can check that the flatness condition Ut — Vx + [U,V] = of this Lax pair yields the 
dAKNS system f lMlISS]) . Note that we can generate a novel two-fold integrable hierarchy for 
this deformed AKNS system. Firstly, by keeping the perturbed equations fl5m52p the same, 
but by increasing the order of the differential constraint recursively as (15^ . we obtain a new 
integrable hierarchy for the dAKNS. Secondlly, by keeping the constraint fixed to its lowest 
level, i.e. as (155|) or (l53ll . one can increase also the order of the AKNS equation with higher 
dispersions in the lefthand side of fl51H52l) in the standard way, generating another integrable 
hierarchy. 

This general deformed AKNS system along with its two-fold integrable hierarchy proposed 
here can yield as particular cases both the dmKdV and the dKdV equations considered above. 
One can check that the dmKdV equation proposed here can be obtained directly from the 
dAKNS through the reduction q = r = v , which degenerates fl5HI52p to the same equation ([3]), 
while the constraint (jl]) is derived from fl55|) with 

G^'^ =tic'it)-w')^a^ + iwa\ G^^^ = 0. (58) 

Similarly for deriving the dKdV from the general case of dAKNS system we have to consider 
the particular reduction q = u,r = 1 , which obviously makes flS^ trivial while reduces fIFI]) 
to the deformed KdV ([T]). At the same time the constraint ([2]) can be derived by solving 

as 

where u{x, t) is the KdV field and (x, t) is the deforming function with c{t) as its asymptotic 
value \im\x\^^ g{x,t) = c(t) . 

Reductions of the deformed AKNS system to other important integrable deformations, e.g. 
NHD of the nonlinear Schrodinger and sine-Gordon equations have been considered in 
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5 Concluding remarks 



We have extended the concept of nonholonomic deformation from the KdV equation to the 
mKdV equation and constructed for this novel deformed system all essential and important 
structures of an integrable system, like a matrix Lax pair, exact N-soliton solution through 
1ST, integrable hierarchies etc. Interestingly, the integrable deformed mKdV, that we have 
discovered here, shows even a richer picture than that of the usual integrability. 

In particular, it allows accelerating exact solitons for the basic field as well as for the 
perturbing function. The perturbing function, similar to the deformed KdV case [3], takes 
a consistent solitonic form and forces the field soliton through its asymptotic value c{t) at 
space infinities, to move with an acceleration (or deceleration). Therefore such a self consistent 
solitonic perturbation, which preserves the integrability, could be significantly important in 
laboratory experiments with mKdV solitons, where the usual loss of energy inevitable in a 
realistic system could be compensated for by the driving force given from the boundaries as 
found here. 

Moreover, the deformed mKdV system shows a novel two-fold integrable hierarchy, as found 
also for the deformed KdV |3] . The first one corresponds to the standard mKdV hierarchy with 
higher dispersions but perturbed with a deforming function subjected to a given nonholonomic 
constraint. The second one is a new hierarchy of equations, where the same mKdV equation is 
perturbed now by a function with increasingly higher order of nonholonomic constraint. 

In parallel with the established procedure for finding the symmetries of the integrable 
equations, we have constructed for the first time an infinite number of generalized symmetries, 
conserved quantities and a recursion operator for both the deformed KdV and the deformed 
mKdV equations. It reveals a remarkable fact that, though the deforming functions are con- 
tained in local current densities, the local densities do not have any dependence on them. Hence 
the integrated conserved quantities are also not influenced by the deformation and remain the 
same as in the unperturbed case. 

Finally we have unified and generalized the deformed KdV and mKdV equations to the 
deformed AKNS system, with explicit construction of its Lax pair and a two-fold integrable hi- 
erarchy. Particular reductions of this new deformed integrable AKNS system yield the deformed 
KdV as well as the deformed mKdV system presented here. 

The study of the deformed Painleve class of equations obtained as reductions of the dKdV 
and dmKdV equations, a potentially important project, we wish to pursue further. 
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FigTirc 1: Exact soliton solutions in the (x,t)-planc: a) for the standard mKdV equation with 
constant soliton velocity (indicated by its constant inclination), b) for the nonholonomically 
deformed dmKdV equation with accelerated soliton (reflected in its curvature with c{t) = Q.5t). 
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